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Abstract. The many-fermion Lagrangian which includes the 't Hooft six-quark flavor mixing interaction 
{Nf — 3) and the Ul (3) x Ur{S) chiral symmetric four-quark Nambu - Jona-Lasinio (NJL) type interactions 
is bosonized by the path integral method. The method of the steepest descents is used to derive the effective 
quark-mesonic Lagrangian with linearized many-fermion vertices. We obtain, additionally to the known 
lowest order stationary phase result of Reinhardt and Alkofer, the next to leading order (NLO) contribution 
arising from quantum fluctuations of auxiliary bosonic fields around their stationary phase trajectories (the 
Gaussian integral contribution). Using the gap equation we construct the effective potential, from which 
the structure of the vacuum can be settled. For some set of parameters the effective potential has several 
extrema, that in the case of SU{2)i x U{1)y flavor symmetry can be understood on topological grounds. 
With increasing strength of the fluctuations the spontaneously broken phase gets unstable and the trivial 
vacuum is restored. The effective potential reveals furthermore the existence of logarithmic singularities 
at certain field expectation values, signalizing caustic regions. 
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1 Introduction 

The global 1/^(3) x Ur{3) chiral symmetry of the QCD La- 
grangian (for massless light quarks) is broken by the Ua{^) 
Adler-Bell-Jackiw anomaly of the SU{3) singlet axial cur- 
rent qjfijsq. Through the study of instantons 021, it has 
been realized that this anomaly has physical effects with 
the result that the theory contains neither a conserved 
U{1) quantum number, nor an extra Goldstone boson. In- 
stead, effective 2Nf quark interactions arise, which are 
known as 't Hooft interactions. In the case of two fla- 
vors they are four-fermion interactions, and the resulting 
low-energy theory resembles the original Nambu - Jona- 
Lasinio model In the case of three flavors they are 
six-fermion interactions which are responsible for the cor- 
rect description of rj and 77' physics, and additionally lead 
to the OZI-violating effects [HtS] . 



chiral symmetric four-quark interaction 



G, 



(2) 



C 



2Nf 



K(det qPLq + det qPuq) 



(1) 



has been additionally included to the quark Lagrangian 

Ant — ^6 + ^4- (3) 

To bosonize the theory in both mentioned cases one has to 
integrate out from the path integral a part of the auxiliary 
degrees of freedom which are inserted into the original 
expression together with constraints [H| 

1 = / Yl'DSa'DpaSiSa ~ q\a(i)^{Pa " qijbKq) 
•' a 

Y[VSaVpaVaaV(t,a (4) 

1 

exp < i / d'^xlaaisa - qXaq) + MPa ^ qilsKq)] 



where the matrices Pi_ij=(l=F 75)72 are projectors and 
the determinant is over flavor indices. 

The physical degrees of freedom of QCD at low-energies 
are mesons. The bosonization of the effective quark inter- 
action (Q) by the path integral approach has been con- 
sidered in . A similar problem has been studied by 
Reinhardt and Alkofer in 8 , where the Ul{3) x Ur{3) 
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The auxiliary bosonic fields, Ua, and, 0^, (a = 0, 1, . . . , 8) 
become the composite scalar and pseudoscalar mesons and 
the auxiliary fields, Sa, and, Pa, must be integrated out. 
The standard way to do this is to use the semiclassical or 
the WKB approximation, i.e. one has to expand the Sa 
and Pa dependent part of the action about the extremal 
trajectory. Both in 6,7 and in 8 the lowest order station- 
ary phase approximation (SPA) has been used to estimate 
the leading contribution from the 't Hooft determinant. In 
this approximation the functional integral is dominated 
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by the stationary trajectories rg-^{x), determined by the 
extremum condition 6S{r) = of the action S{r)^. The 
lowest order SPA corresponds to the case in which the 
integrals associated with S'^S{r) for the path rg^(x) are 
neglected and only S{rg^) contributes to the generating 
functional. 

In this paper we obtain the ?i-correction to the leading 
order SPA result. It contains not only an extended version 
of our calculations which have been published recently [S] 
but also includes new material with a detailed discussion 
of analytic solutions of the stationary phase equations, cal- 
culations of the effective potential to one-loop order, solu- 
tions of the gap equations, general expressions for quark 
mass corrections, and quark condensates. We also discuss 
the results of the perturbative approach to find solutions 
of the stationary phase equations. 

There are several reasons for performing the present 
calculation. First, although the formal part of the problem 
considered here is well known, being a standard one-loop 
approximation, these calculations have never been done 
before. The reason might be the difficulties created by 
the cumbersome structure of expressions due to the chiral 
group. Special care must be taken in the way calculations 
are performed to preserve the symmetry properties of the 
theory. Second, it provides a nice explicit example of how 
the bosonization program is carried out in the case with 
many-fermion vertices. Third, since the whole calculation 
can be done analytically, the results allow us to examine 
in detail the chiral symmetry breaking effects at the semi- 
classical level. By including the fluctuations around the 
classical path related with the 't Hooft six-quark determi- 
nant, our calculations of the gap equations and effective 
potential fill up a gap existing in the literature. Fourth, 
the problem considered here is a necessary part of the 
work directed to the systematic study of quantum effects 
in the extended Nambu - Jona-Lasinio models with the 't 
Hooft interaction. It has been realized recently that quan- 
tum corrections induced by mesonic fluctuations can be 
very important for the dynamical chiral symmetry break- 
ing JUJJm although they are 1/iVc supressed. 

Let us discuss shortly the main steps of the bosoniza- 
tion which we are going to do in the following sections. As 
an example of the subsequent formalism, we consider the 
bosonization procedure for the first term of Lagrangian 
Using identity Q one has 



exp — j d^x Kdet{qPLq) 



= / YlT^Waexpf^ / d''x{-qWPLq) 



(5) 



n -DUl exp\^ I d^x ( i W^„C/t + ^ det 



The variables W = WqAq, where Wa — (^a — 'i-4'a, describe 
a nonet of meson fields of the bosonized theory. The aux- 
iliary variables W = UlXa, Ul = Sa + ipa must be inte- 

^ Here r is a general notation for the variables {sa,Pa) and 
S{r) is the r-dependent part of the total action. 



grated out. The Lagrangian in the first path integral as 
well as the Lagrangian L{W, W) in the second one have 
order iV°, because qq and count as Nc, k ^ and 
W ~ Thus we cannot use large Nc arguments to ap- 

ply the SP method for evaluation of the integral. However 
the SPA is justified in the framework of the semiclassi- 
cal approach. In this case the quantum corrections are 
suppressed by corresponding powers of h. The stationary 
phase trajectories are given by the equations 



dL Wa 3k 



dUl 2 64 



AabcUM^O 



(6) 



where the totally symmetric constants, Aabc, come from 
the definition of the flavor determinant: 



det = AabcUlulu^, , 



and equal to 



1 



-^ahc — 21 ^ijk^mnl (-^a )'tm ('^6 )jn ('^c)fc/ 



(7) 



2 / 2 

= -^dabc + Y 3 {SSaoSboScO — ^aoSbc — SbO^ac — ScoSab) 

with Xa being the standard U{3) Gell-Mann matrices, 
[Xa,Xb] = 2ifabcXc, {Xa,Xb} = 2dabcXc, normalized such 
that trAaAfc — 25ab, and a = 0, 1, . . . , 8. 

The solution to Eq.©, U'^^(yV), is a function of the 
3x3 matrix W 



Ul^{W) = 4,\l—W~'^{deiW)^. 



(8) 



Expanding L{W^U^) about the stationary point U}^j^{W) 



we obtain 



L{W, U^) = L{W, uUw)) + ]-ul . ill 
1 d^L 



3! "dUldUldUl ' ' 
where {7^ = [/^ — U^^ and, as one can easily get, 

m, Ul^iW)) = ItriWUl,) + ^ det C/J^ 
det IF 



(9) 



Therefore, we can present Eq.Q in the form 
exp|— J d'^x K det {qPj^q) 

= jYlvWaexpi^^J d^x (^^qWPLq + 2 



(10) 



(11) 



det IF 
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which sphts up the object of our studies in two contri- 
butions which we can clearly identify: the first line con- 
tains the known tree-level result jT, and the second line 
accounts for the ?i-suppressed corrections to it, which we 
are going to consider. Unfortunately, the last functional 
integral is not well defined. To avoid the problem, we 
will study the theory with Lagrangian In this case 
the functional integral with quantum corrections can be 
consistently defined in some region where the field- 
independent part, Dab, of the matrix d^L{r)/dradri, has 
real and positive eigenvalues. In order to estimate the ef- 
fect of the new contribution on the vacuum state we derive 
the modified gap equation and, subsequently, integrate it, 
to obtain the effective potential V{J-). On the boundary, 
dJ^, the matrix Dab has one or more zero eigenvalues, 
(io(/c) = 0, and hence Dab is noninvertible. As a con- 
sequence, the effective potential blows up on dJ^. This 
calls for a more thorough study of the effective potential 
in the neighbourhood of dJ-', since the WKB approxima- 
tion obviously fails here (region of the caustic). Sometimes 
one can cure this problem going into higher orders of the 
loop expansion |12II13| . Nevertheless, V{T) can be ana- 
lytically continued for arguments exterior to dJ-, where 
(io(fc) are negative. In fact, because of chiral symmetry, we 

have two independent matrices Z?^),'' and D^^J associated 

with the two quadratic forms SaD^Jjj Sb and paD^^Jpb in 
the exponent of the Gaussian integral. The eigenvalues 
of these matrices are positive in the regions J-i and J-2 
correspondingly, and !Fi D ^2- Accordingly, the SU{3) 
effective potential is well defined on the three regions: 
= {d^o\So^ > 0}, ^2 - {4^' > 0,d^o^ < 0}, and 

•^3 = {'^o^^'^o^'' < 0} separated by two boundaries 
and dTi where V +oo. It means that the effective 
potential has one stable local minimum in each of these 
regions. However, we cannot say at the moment how much 
this picture might be modified by going beyond the Gaus- 
sian approximation near caustics. 

Our paper is organized as follows: in Sec. |2 we de- 
scribe the bosonization procedure by the path integral for 
the model with Lagrangian (O and obtain h corrections to 
the corresponding effective action taking into account the 
quantum effects of auxiliary fields Va ■ We represent the La- 
grangian as a series in increasing powers of mesonic fields, 
aa,(t>a- The coefficients of the series depend on the model 
parameters G, k, m, and are calculated in the phase where 
chiral symmetry is spontaneously broken and quarks get 
heavy constituent masses m„,TOd,m.s. We show that all 
coefficients are defined recurrently through the first one, 
ha- Close- form expressions for them are obtained in Sec. 13 
for the equal quark mass as well as m„ — nid ^ cases. 
In Sec. Ol we also study U corrections to the gap equa- 
tion. We obtain ?i-order contributions to the tree-level con- 
stituent quark masses. The effective potentials with SU{3) 
and SU{2)j x U{1)y flavor symmetries are explicitly cal- 
culated. In Sec. 01 we alternatively use the perturbative 
method (1/iVc-expansion) to solve the stationary phase 
equations. We show that this approach leads to strong 
suppression of quantum effects. The result is suppresed 



by two orders of the expansion parameter. We give some 
concluding remarks in Sec. El Some details of our calcula- 
tions one can find in three Appendices. 



2 Path integral bosonization of many-fermion 
vertices 

The many-fermion vertices can be linearized by introduc- 
ing the functional unity Q in the path integral represen- 
tation for the vacuum persistence amplitude 



Z = VqVqexp 



(12) 



We consider the theory of quark fields in four dimen- 
sional Minkowski space, with dynamics described by the 
Lagrangian density 



C = q{ij''df, - m)q + £5^^. 



(13) 



We assume that the quark fields have color {Nc — 3) 
and flavor (Nf = 3) indices which range over the set 
i = 1,2,3. The current quark mass, m, is a diagonal matrix 
with elements diag(m„, to^, jtis), which explicitly breaks 
the global chiral S'?7l(3) x SUu{3) symmetry of the La- 
grangian. The second term in (|13(l is given by (j^J. 

By means of the simple trick (0J), it is easy to write 
down the amplitude (fT^ as 

„ 8 8 8 8 

Z= VqVq Jl VSa \{ Vpa \{ VcJa n '^(^a 
a=a a=0 a=0 a=0 

exp (i f d'^xC] (14) 



with 



G 



C' = q{i-f^d^ ~m-a- i-fr,^)q + - [(sa)' + (Pa)^] 
+ Sa(7a + Pa4>a + ^ [det(s + ip) + det(s - ip)] (15) 

where, as everywhere in this paper, we assume that a — 
(Ja^a, and so on for all auxiliary fields: (/>, s, p. Ea. H14|l 
defines the same expression as Ea. (|12f) . To see this, one 
has to integrate first over auxiliary fields <7a, (j)a- It leads 
to (5-functionals which can be integrated out by taking in- 
tegrals over Sa and Pa and which bring us back to the 
expression (I12f) . From the other side, it is easy to rewrite 
Ea. (|14(l . by changing the order of integrations, in a form 
appropriate to accomplish the bosonization, i.e., to calcu- 
late the integrals over quark fields and integrate out from 
Z the unphysical part associated with the auxiliary Sq, pa 
bosonic fields, 



Z = J Y['D<Ta'D(t)a'Dq'Dqexp J d'^xCq{q,q,<J,( 

J Y\_'^Sa'Dpaexp (^i J d'^xCr{a,4>,s,p)^ (16) 
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where 

l^q = q{il^df^ - m - cr - ij5(t>)q, (17) 
G 

= [{Sa)^ + (Pa)^] + (SaCa + Pa4>a) 

+ -^AabcSa (SfcSc " iPbPc) ■ (18) 

The Fermi fields enter the action bilinearly, we can always 
integrate over them, because in this case we deal with 
a Gaussian integral. At this stage one should also shift 
the scalar fields aa{x) — > cra{x) + Aa by demanding that 
the vacuum expectation values of the shifted fields vanish 
< 0|(Ta(a;)|0 >= 0. In other words, all tadpole graphs in 
the end should sum to zero, giving us the gap equation to 
fix constants Aa- Here Aa = ma — rha, with nia denoting 
the constituent quark masses ^ . 

To evaluate functional integrals over Sa and Pa 

+ OC 



/ + 00 
YlVSaVpa 
-oo 



expl^i J d^xCria + A,<j),s,p)j (19) 

where Af is chosen so that Z[0, 0; A] ~ 1, one has to use 
the method of stationary phase. Following the standard 
procedure of the method we expand Lagrangian Cr{s,p) 
about the stationary point of the system r^^ = (s^^, Pg^)- 
Near this point the Lagrangian Cr{s,p) can be approxi- 
mated by the sum of two terms 

£^((7 + A,(l),s,p) sa Cr{rst) 

+ ^E^~"(^)'^"/3(^st)^>(^) (20) 

a, (3 

where we have only neglected contributions from the third 
order derivatives of Cr{s,p). The stationary point, r^^, is 
a solution of the equations C'j.{s,p) = determining a flat 
spot of the surface £r{s,p): 



3k 

GSa + (o- + A)a + —Aabc{sbSc " PbPc) = 



(21) 



Gpa + (l^a- ^AabcSbPc 

16 



0. 



This system is well-known from We use in Ea. (|20|l 
symbols f" for the differences (r" — r^^). To deal with the 
multitude of integrals we define a column r with eighteen 
components fa = {sa > Pa) and with the real and symmetric 
matrix C^^irgi) being equal to 



^a/3(^St) 



GSab + ^AabcSl^ -^AabcPlt 
-.AabcP^l GOab ~ —AabcS. 



16 



16 



St 



(22) 

^ The shift by the current quark mass is needed to hit the 
correct vacuum state, see e.g. |14| . The ftmctional integration 
measure in Eg. 11611 does not change under this redefinition of 
the field variable aa{x). 



The path integral p9(l can now be concisely written as 
Z[a,(/);Z\] « AAexp (^i J d^xCAr^t)^ J^W^ 

expQ j d'xf\x)L';{r^^)f{x)y (23) 

The Gaussian multiple integrals in Ea. (|23|l define a func- 
tion of Ca/jir^i) which can be calculated by a general- 
ization of the well-known formula for a one-dimensional 
Gaussian integral. Before we do this, though, some addi- 
tional comments should be made: 

(1) The first exponential factor in Ea. l|23|l is not new. 
It has been obtained by Reinhardt and Alkofer in 0. A 
bit of manipulation with expressions H18|l and H21|l leads 
us to the result 

^^(^St) = \ {G[{slO^ + [Pltf] + 4[(a + A)aSl^ + <j^aPlt]} 



- Y^tr(f/stC/s\) + t^^^Kl + ^'U, 



St) 



(24) 



Here the trace is taken over fiavor indices. We also use the 
notation W — WaXa and U — UaXa where Wa = aa+Aa — 
iff'a, Ua — Sa — ipa- H is similar to the notation chosen in 
Eq.© with the only difference that the scalar field aa is 
already splitted as (Ja o'a + Aa- This result is consistent 
with ifTUI) in the hmit G — 0. For this partial case Ea. lf^ 
coincides with Eq. © and we know its solution JHJ . If G 7^ 
we have to obtain the stationary point Ug^ from Eq. (|21|) . 

One can try to solve Eqs.|2U exactly, looking for solu- 
tions Sg|. and in the form of increasing powers in fields 



Sgt =ha + /lifcVfc + /libcCTbO-c + h^abcMc 

+ h^abcd(^b<Tc(Td + h''abcd'^b(l)c(t>d + ■■■ (25) 

Pst = '^aUb + h^abc^t'b^c + h^abcd'^b(Tc(t>d 

+ h^abcdMc<Pd + . . . (26) 

with coefficients depending on Aa and coupling constants. 
Putting these expansions in Eqs.|2U obtains a series 

of selfconsistent equations to determine ha, h^^J , /i^^"* and 
so on. The first three of them are 

3k 



Gha + Aa + —Aabchhc = 0, 



Gdac 
GSac 



3k 
l6 
3k 

Te' 



Anrhhb I hi 



(27) 



Aacbhb I h^c} 



All the other equations can be written in terms of the 
already known coefficients, for instance, we have 



'abc 



' (3) 
abc 



3k 



hnn h,T hr-r- A 



32 bb 



cc ^abc ' 



3k 
16' 



;,(2), (2), (1) . _ 
"aa "'bb cc -^abc i 



^^^l^-^-h^lh^^^h^^A-ab,, 
aoc oo Q", cc aoc ' 

„3k (1) (1) (1) 



abed 



-|^g "aa "bb "'ccd^abc 
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"•abed 



(1) fuWi,{2) 7,(2);, (3) \ A 



It is assumed that coupling constants G and k are chosen 
such that Eas. l|27|l can be solved. Let us also give the rela- 
tions following from H27() which have been used to obtain 



{GK + 2zi„)/iiV = -{iGha + 2Aa)h'- 



(2) 
b ■ 



(29) 



As a result the effective Lagrangian (|^ can be expanded 
in powers of meson fields. Such an expansion, up to and 
including the terms which are cubic in aa,4>a, looks like 



,(2), 



©(field"*). 



(30) 



This part of the Lagrangian is responsible for the dynam- 
ical symmetry breaking in the quark system and for the 
masses of mesons in the broken vacuum. 

(2) The coefficients ha are determined by couplings 
G, K and the mean field Z\q. This field has in general only 
three non-zero components with indices a = 0, 3, 8, ac- 
cording to the symmetry breaking pattern. The same is 
true for ha because of the first equation in (|27|l . It means 
that there is a system of only three equations to determine 
h = haXa = diag(/i„, hd, hs) 



Ghu + Au + -Tihdhs = 
16 



Ghd + Ad + -^hshu = 
16 



Ghs + As + —huhd = 
16 



(31) 



This leads to a fifth order equation for a one-type variable 
and can be solved numerically. For two particular cases, 
rhu = rhd = rhg and rha — rhd ^ w^, Eas. (|31|l can be 
solved analytically, because they are of second and third 
order, correspondingly. We shall discuss this in the next 
section. 

(3) Let us note that Lagrangian (|18|l is a quadratic 
polynomial in Pa{x) and cubic with respect to Sa{x). It 
suggests to complete first the Gaussian integration over 
Pa and only then to use the stationary phase method to 
integrate over s^. In this case, however, one breaks chiral 
symmetry. This is circumvented by working with the col- 
umn variable Tq,, treating the chiral partners {sa,Pa) on 
the same footing. It is easy to check in the end that the 
obtained result is in agreement with chiral symmetry. 

Let us turn now to the evaluation of the path integral 
in Ea. (|23(l . After the formal analytic continuation in the 



time coordinate xq — s- 1x4, we have'^ 

/+00 
-CO 



exp 



i J d''xEr\xE)C':ir^t)r{xE?j , (32) 



where the subscripts E denote Euclidean quantities. To 
find an expression for -fi'[''st]; split the matrix £" into 
two parts £"(rg^) = D — L where 



h 



a/3 



h 





(2)-l 
ab 



3k 

^al3 — ~r^Aabc 

16 



st ■ 
Pit 



/ Q/3 

he) p. 
*st 



St 



(33) 



Q/3 



The matrix D corresponds to evaluated at the point 
ha, Pqf — and is simplified with the help of 



'st 



Eas. H27|l . The field-dependent exponent with matrix L can 
be represented as a series. Therefore, we obtain 

^['^St] ^ j W'DfaGX^ j dl^XEfaDaprp 

X];;^ ( 2 / ^'^^^^"^"''^'3 ) • (34) 



The real symmetric matrix D^p can be diagonalized by 
a similarity transformation = Sapr'^. We are then 
left with the eigenvalues of the matrix D in the path 
integral H34I) . These eigenvalues are real and positive in 
a finite region fixed by the coupling constants. For in- 
stance, if rhu — "hid = "his the region is (— kZ\„) < 12G^ 
where, as usual in this paper, we assume that G > and 
K < 0. In this region the path integral is a Gaussian 
one and converges. To perform this integration, we first 
change variables = Upaqcr, such that U rescales eigen- 
values to 1. The quadratic form in the exponent becomes 
faDapfp = f'^{S* DS)ai3f'i^ — (lala- The matrix of the to- 
tal transformation, Vaa — SapUpa, has the block-diagonal 
form 



/F^'^ 



"-ah ^ *^ac ^bc ' 



ap 



(35) 



Then the integral (|34|l can be written as 

K[rst] = <^^ctV J JIX'Qaexp ^-^ J d-^XEQaqa 
00 . . „ 

J2^i^-jd''xEqc.{V'LV)e,pqpj . (36) 



It differs from the standard Wick rotation by a sign. The 
sign is usually fixed by the requirement that the resulting Eu- 
clidean functional integral is well defined. Our choice has been 
made in accordance with the convergence properties of the path 
integral H23|l . 
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By replacing the continuum of spacetime positions with 
a discrete lattice of points surrounded by separate re- 
gions of very small spacetime volume fi, the path inte- 
gral H36(l may be reexpressed as a Gaussian multiple in- 
tegral over a, finite number of real Vciriables Qa.,xi 

where 

-°° a,x \ a,x J 



n=0 



(37) 



where the matrix W is given by yVax,i3y — 

The Gaussian integrals in this expression are well known 



Yl dqt [qki 9/02 ■• ■ 9fc2„ ) exp [ - y ^ Qi 



1 



(38) 



Here Sk^k^.-Mn ^ totally symmetric symbol which gen- 
eralizes an ordinary Kronecker delta symbol, 5ij , with the 
recurrent relation 

Skik2...k2„ — '5felfe2'^fcife2...fc2„ + '^fclfe3'^fcife2fe3...fc2„ 

+ Sk^k2Jk,k2..M,.■ (39) 

The hat in this formula means that the corresponding 
index must be omited in the symbol S. Let us also re- 
mind that integrals of this sort with an odd number of 
g-factors in the integrand obviously vanish. The multiple 
index k is understood as a pair k — a.x with the Kro- 
necker 4ifc2 = ^0102^X1X2- 

By performing the Gaussian integrations one can fi- 
nally fix the constant of proportionality J\f and find that 

1 



^[^St] = 1+ 2'5fclfe2Wfeife2 + ■•• 

+ ^^<5fclfc2...fc2„ Wfcifc2 . . . Wk2 



lk2n 



.(40) 



The infinite sum here is nothing else than 1 / y'Det (1 — W). 
The determinant may be reexpressed as a contribution to 
the effective Lagrangian using the relation Det (1 — W) = 
expTrln(l — W). Thus, we have 



K[rst] =exp --Trln(l-W) 



(41) 



with the logarithm of a matrix defined by its power series 
expansion 

ln(l - W) = -W - - - ... 



b^^y (ln(l - V^LY)) 



a,f3 ■ 



(42) 



The path integral (|34|l defines a function of Dap that is 
analytic in -Da/3 in a region around the surface where the 



eigenvalues of Dap are real positive and the integral con- 
verges. Since Ea. H41|) equals to (|34|) . it provides the ana- 
lytic continuation of Ea. H34l) to the whole complex plane, 
with a cut required by the logarithm. 

Let us now return back from the spacetime discret 
lattice to the spacetime continuum. For that one must 
take the limit — > oo in Ea. H41|l . and replace . . . — + 
J d*XE .... As a result we have 

K[rst] = exp (^-^^"' J d^^E tr ln(l - V'LV)^ 

= exp J d'^XE tr ln(l + F)^ , (43) 

where 'tr' is to be understood as the trace in an ordinary 
matrix sense. In the second equality we have used the 
property of matrix V, given by Ea. H35|) . This property 
can be used because of the trace before the logarithm. 
The matrix F is equal to 



3k 

i'af} = -—Ache 

16 



Q/3 



The factor fl ^ may be written as an ultraviolet divergent 
integral 12"^ = 5%{x — x). This singular term needs to 
be regularized, for instance, by introducing a cutoff Ae 
damping the contributions from the large momenta kE 



n-^ = 4(0) ^ ^ 



E A% 



^,/2 (2^)4 (2^)4 



(45) 



To finish our calculation one needs to return back to 
the Minkowski space by replacing —ixQ. It follows 

then that the functional integral (|23|l is given by* 



Z[a,0;Z\] 



(46) 



The action (|46|l contains in closed form all informa- 
tion about 7i-order corrections to the classical Lagrangian 
Cr{rgf^). Nevertheless it is still necessary to do some work 
to prepare this result for applications. In the following we 
deal mainly with the first term of the series, since it is the 
only one that contributes to the gap equation. 



St 



(47) 

Here one should sum over indices a,6, c ~ 0,1,... 8. In 
the next section we will calculate these sums for the cases 
with exact 51/(3) and broken S'C/(3) ^ SU{2)i x U{1)y 
flavor symmetry. 

To give some additional insight into the origin of for- 
mula let us note that ?i-order corrections to £r(''st) 

* The sign of the term ~ must be corrected accordingly 
in Eq.(34) of 0. 
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can be obtained without evaluation of the Gaussian path 
integral in Eq. H23() . Instead one can start directly from La- 
grangian and obt ain the one-loop contribution using 
the canonical operator formalism of quantum field theory 
|15) . Since the Lagrangian does not contain kinetic terms 
the time-ordered products of meson fields have a pure sin- 
gular form: < T{aa{x)ab{y)} >= iS{x - yWab~^ ' < 

T{(j)a{x)(j)b{y)} >= iS{x - y)h\^ . It is easy to see, for 
instance, that the tadpole contribution, coming from the 
cubic terms in (|3()|l . exactly coincides with the term linear 
in a in Ea. lHTjl . 

3 The ground state in the semiclassical 
expansion 

The considered model belongs to the NJL-type models 
and therefore at some values of coupling constants the 
many-fermion interactions can rearrange the vacuum into 
a chirally asymmetric phase, with mesons being the bound 
states of quark and antiquark pairs and unconfined quarks 
with reasonably large effective masses. The process of the 
phase transition is governed by the gap equation, and, as 
we already know from Eq. H47|l , the gap equation is modi- 
fied by the additional contribution which comes from the 
term ~ f2^^ . Our aim now is to trace the consequences of 
this contribution on the ground state. 

An effective potential V{aa) that describes the system 
in the chirally asymmetric phase has a minimum at some 
non-zero value of ctq = Aa- The effective Lagrangian con- 
structed at the bottom of this well does not contain linear 
terms in a fields. It means that the linear terms in Ea. (|47|l 
must be canceled by the quark tadpole contribution. This 
requirement can be expressed in the following equation 

hi + ■^miJo{mf) = Zi , 

(48) 

The relation of flavor indicies i = u,d,s and indicies a, 6, c... 
is given below in Ea. (|51|l . The first term, hi, and the term 
on the right-hand side, Zi, are the contributions from La- 
grangian H47|l . The second one is the contribution of the 
quark loop from H17() with a regularized quadratically di- 
vergent integral Jo(to^) being defined as 

= -m^\n(l + ^\ . (49) 

This integral is a positive definite function for all real val- 
ues of the cutoff parameter A and the mass m. The kernel 
p{t, yl^) = 1 - (1 -I- tA^) exp(-tyl2) is introduced through 
the Pauli - Villars regularization of the integral over t, 
which otherwise would be divergent at the point t = 0. 
Assuming that hi is a known solution of Eas. (|31|l . being a 
function of parameters G, k and the vacuum expectation 



values Ai of scalar fields, we call Ea. H48|l a gap equation. 
The right-hand side of this equation is suppressed by a 
factor h in comparison with the first two terms. Thus an 
exact solution to this equation, i.e. the constituent quark 
masses rrii {i = u, d, s), would involve all powers of h, but 
higher powers of h in such a solution would be affected by 
order h^ corrections to the equation. It is apparent that 
one has to restrict our solution only to the first two terms 
in an expansion of rrii in powers of h, obtaining rm in 
the form rrii = Mi + Ami. It is important to stress for 
physical applications that Ea. (|48|) gets at h order an ad- 
ditional contribution from meson loops |l()llllj which we 
do not consider here. We restrict our attention only to the 
new kind of contribution of order h which is essential for 
the Nf = 3 case and has not yet been discussed in the 
literature. 



3.1 Gap equation at leading order: SU{3) case, two 
minima 

The first two terms of Eq. H48|l are of importance at lead- 
ing order. Combining this approximation for Ea. l|48() (i.e. 
setting Zi = 0) together with Ea. (pTT|l . one obtains the 
gap equation already known from a mean field approach 
|3f IB' , which self-consistently determines the constituent 
quark masses Mi as functions of the current quark masses 
and coupling constants 

Ai — —Ghi — — tijkhjhk, 

hi = --:r-^miJo{mf) = 2 < q^i >= 2ai{mi). (50) 

Here the totally symmetric coefficients tijk are equal to 
zero except for the case with different values of indices 
* J 7^ ^ when tuds = 1- The latin indices i,j,k mark 
the flavour states i = u,d, s which are linear combinations 
of states with indices 0,3 and 8. One projects one set to 
the other by the use of the following matrices ujia and Cai 
defined as 

= ' = ^ 1 

1 1-2/ \ V2 -2/ 

(51) 

Here the index a runs a = 0,3,8 (for the other values of 
a the corresponding matrix elements are assumed to be 
zero). We have then ha = Saihi, and hi = ujiaha- Similar 
relations can be obtained for Ai and Aa- In accordance 
with this notation we use, for instance, that = ujiah^ca ■ 
The following properties of matrices (|51|l are straightfor- 
ward: UOiaeaj = Si-j, CaiUOib = Sab and CaiCaj = %/2. The 

coefficients i^jfe are related to the coefficients Aabc by the 
embedding formula ^lOJiaAabcGbje-ck = Ujk- The SU{3) ma- 
trices Aq with index i are defined in a slightly different way 
2Xi = LOiaXa and — 2eaiXi. In this case it follows that, 
for instance, a = (Ja^a = CTjAi ~ diag((7„, dd, CTs), but 

2aaAa = (TiAi. 
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Fig. 1. The function hu{mu) (here [hu] = GeV^, [niu] = GeV) 
at fixed values m„ = 6 MeV, A = 860 MeV for G = 0, k = 
-5000 GeV"'^ (large dashes), G = 5 GeV"^ k = (dash 
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Fig. 2. The h leading order effective potential V{mu) corre- 



dotted line), G = 5 GeV-^^ = -5000 GeV'^ (solid line). sponding to the pattern s^hown by the solid line on FigE] We 
The small-dashes hue corresponds to the function 2a„(m„). ^^'^ ^^^^^ t'^l - I"^"] ~ 



It is well known 0] that the gap equation at leading 
order has at least one nontrivial solution, Mi{'mj, G, k, A), 
if K, rhi 7^ 0. This important result is a direct conse- 
quence of the asymptotic behavior of the quark condensate 
ai{mi) : ai{0) = ai{oo) — 0, ai{mi) < if < < cxd, 
from one side, and the monotonical decrease of the neg- 
ative branch of the function hi{Aj) on the semi- infinite 
interval rhi ^ "Tfii < cxd, from the other side. Stronger 
statements become possible if we have more information. 
Let us assume that SU(3) flavour symmetry is preserved. 
If rhu = rhd — rhs, one can conclude that Z\„ = Ad — 
and, as a consequence, we have hu = hd~ hg ~ ^2/3/io- 
Instead of a system of three equations in H5U|) there is only 
one quadratic equation with two solutions 




(52) 



The second solution is always positive (for k < 0, G > 0) 
and, due to this fact, cannot fulfill the second equation in 
(jSUIl . On the contrary, the solution hu \ as k ^ 0, gives 
hu^ = —Au/G, which leads to the standard gap equation 
2ir^Ai = NcGmiJo{mi) for the theory without the 't 
Hooft determinant. Alternatively one could consider the 
theory with only the 't Hooft interaction Q taking the 

limit G ^ in Ea. (|52|l . In this case we obtain for hi}'' 



hi^'> = -V-kAu = 4 sign(K) J ^ 

K V 



(G = 0). (53) 



In Fig^we plot three different stationary trajectories 
hu — hu(jnu) (with corresponding parameter values given 
in the caption) as functions of the quark mass to„. To be 
definite we put the current quark mass equal to„ = 6 MeV 
and the cutoff parameter is choosen to be yl = 860 MeV. 
By fixing A we completely fix the curve corresponding to 
the right-hand side of the gap equation, i.e. the function 
2q!„ = —NcrauJoira^) I {2i{^). This function starts at the 
origin of the coordinate system, being always negative for 
positive values of m^. At the point m„ = m ~ 0.74/1 it 
has a minimum ai"""' = —A'^fh/[i{A^ + 'm?')\. All station- 
ary trajectories cross the TO„-axis at the point m„ = rhu- 
The straight line corresponds to the case n — Q. The sec- 
ond limiting case, G = 0, is represented by the solution 
(|53ll and marked by G = 0. The solid curve corresponds to 
Eq. (|52(l , starts at the turning point A with the coordinates 
A{mu,hu) = {rhu + 4G^/k, — 8G/k) and goes monotoni- 
cally down {hu'') and up (/i^f^) for increasing values of 
rriu- The standard assignment of signs for the couplings G 
and k: G > 0, At < is assumed. The points where /iu(m„) 
intersects 2q;„ are the solutions of the gap equation. Let us 
remind 0] that there are two qualitatively distinct classes 
of solutions. The first one is known as a solution with 
barely broken symmetry. We have this solution when /i„ 
intersects 2a„ on the left side of its minimum (the mini- 
mum of a„ is outside the region shown in the figure). If 
hu crosses 2au on the right side of its minimum it cor- 
responds to the case with the firmly broken symmetry. 
It is known that the theory responds quite differently to 
the introduction of a bare quark mass for these two cases. 
The barely broken regime is characterized by strong non- 
linearities reflected in the behaviour of expectation values 
of the scalar quark densities, qiQi, in the physical quark 
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states. Nevertheless, the solutions with the barely broken 
symmetry are likely to be more reliable from the physi- 
cal point of view, in particular, when A ~ 0.5 GeV and 
higher. 

The six-quark interactions add several important new 
features into the picture. For instance, for some set of pa- 
rameters, when 4G^ ~ JtImIkI or G = 0, one can get three 
solutions of the gap equation, instead of one when k = 0. 
One of these cases is illustrated in Fig^Jfor the parameter 
set G = 5 GeV"^ k = -5000 GeV"^ The first solution 
is located quite close to the current quark mass value and, 
being a minimum of the effective potential, corresponds 
to the regime without the spontaneous breakdown of chi- 
ral symmetry. The next solution is a local maximum. The 
third one is a minimum and belongs to the regime with 
barely broken phase. The types of extrcma arc shown in 
FigEl where we used Ea. (|7?1| with J7^^ = for the effec- 
tive potential at leading order. 

The second new feature is that the stationary trajec- 
tories huirUu) in the case k 7^ are real only starting from 
some value of m„ > rrimin- For other values of where 
rriu < rrimin the effective energy is complex. To exclude 
this unphysical region the effective potential V{mu) must 
be defined as a single- valued function on an half-open in- 
terval m„ > m^in. 



3.2 NLO corrections: Zi in the case of SU{3) symmetry 

The functions Zi can be evaluated explicitly. To this end 
it is better to start from the case with SU (3) flavor sym- 
metry. This assumption involves a significant simplifica- 
tion in the structure of the semiclassical corrections, giv- 
ing us however a result which possesses all the essential 
features of the more elaborate cases. Noticing that as a 
consequence of SU{3) symmetry the function ha has only 
one non zero component ho, one can find 

Aacbhb = AacohQ = - {25aO&cO - SaaSca) hu (54) 

where a = 1, 2 ... 8. This result determines in general the 
structure of all flavor multi-indices objects like, for in- 
stance, h'^ab'^^ defined by the Ea. l(?7jl . Indeed, taking into 
account H54(l . one can represent algebraic equations for 

^Ih'^^ in the form 



G [SaoScO (1 ± 2w) -I- SaaSca (1 T t^)] K 



(1,2) 
cb 



^ab , 



16G 



(55) 



where we follow the notation explained in Appendix ^ 
Since this expression is a diagonal matrix, one can easily 
obtain 



Jl,2) 
''cb 



ScoSl 



cOObO 



ScaSi 



ba 



G(l±2w) G(lTw) ' 



(56) 



with associated with the upper sign. 



Consider now the perturbation term of Ea. H48|l . One 
can verify that 



/ 



("ab - n-ab J^abc 



8CP 



2A 



OOc Q=l 



1 - 4w2 1 



V 



-ScQ- 



2G2 V 3 '"(l-4a)2)(i„^2) ■ 
Here we used the properties of the coefficients Aabc 

8 



2 2 



[2 _ ' [2 

\ o ^cO, ^abb = ^abb — / -SaQ 

^ b=0 ^ 



(57) 



(58) 



The quantum effect of auxiliary fields to the gap equation 
stems from the term 



ih['^-h'-:^Mabch 



(1) 



l-3c^2 



3G3 (l + 2w)(l-4w2)(i_^2) 

(59) 

where the right-hand side is the same for the three possible 
choices of the index i — u,d, s. We conclude that in the 
SU{3) limit the functions Zi are uniquely determined by 



UJ{1 - 3^2) 



2G2 {1 + 2lu){1 ~ 4lu^){1 ~ u^) 



(60) 



A direct application of H60|) is the evaluation of the 
corrections Ami to the constituent quark mass. 



^^- = ^^ (1-4C.2)(1_..2) 2^ 



(61) 



This result is based on Eq. (|121l) and Appendix ^ The 
right-hand side must be calculated with the leading order 
value of m„ = M^- The factor Q„ is 



1 



l-(l + 2w)/„(m2) 



(62) 



Here the function Ii{rnf) is proportional to the derivative 
of the quark condensate < qiQi > with respect to 

N G 

h = -^[Jo{m'^)-2m'^Ji{mf)], i = u,d,s (63) 
where the integral Ji is given by 



p{t,A^) 



In 1 



A^ 



A^ 



(64) 



One can show that Qu is related to the expectation val- 
ues of the scalar quark density, qiqi in the physical quark 
state \ Qu>- To be precise we have Q„ = Quu + Qud + Qus, 
where the expectation values Qui are given in Appendix 
IbI It is of interest to know the sign of the quasi-classical 
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(65) 



correction Amu- In general the answer on this question de- 
pends on the values of coupling constants which should be 
fixed from the hadron mass spectrum. All what we know 
at the moment is only that G > 0, k < 0, rhu — 6 MeV 
and yl ~ 1 GeV. One can expect also that the dynami- 
cal masses of the quarks M„ are close to their empirical 
value Mu ^ Mn/3 ^ 300 MeV, with Mn the nucleon 
mass. Let us suppose now that < w <C 1, what actually 
means that the coupling constants belong to the interval 
< —k{Mu — m„)/(4G)^ ^ 1. This range is preferable 
from the point of view of l/Nc counting. In this case we 

" ~ 2G 

and one can conclude that the sign of Aruu is opposite to 
the sign of Q%^^ ■ In turn the function Q^^° is positive in 
some physically preferable range of values of G such that 
G~ 5 GeV-^ 

It must be emphasized that the approximation made in 
Eq.jnU is legitimate only if the quasi-classical correction 
Arriu is small compared with the leading order result Af„. 
In particular, it is clearly inapplicable near points where 
the function Quirriu) has a pole. One sees from Ea. H62|l 
that this takes place beyond some large values of G or 
K. There is a set of parameters for which the function 
(1-1- 2uj)Iu{m'^) is always less than 1. This is the case, for 
instance, for the choice just considered above. For large 
couplings, Qu may have a pole, and one has to check that 
the mean field result M„ is located at a safe distance from 
them before using formula H61() . The large couplings con- 
tain also the potential danger to meet the poles at the 
points uj = 1 and uj = 1/2. These poles are induced by 
caustics in the Gaussian path integral and occur as singu- 
larities in the effective potential. 

We conclude this section with the expression for the 
quark condensate at next to leading order in h, 

< quQu >=< quQu >o 

where the subscript denotes that the expectation value 
has been obtained in the mean field approximation. 
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Fig. 3. The effective potential l/(m„) to ft-order corresponding 
to the case k = -1800 GeV"^ G = 10 GeV"^ rh = and 
A = 860 MeV, where the values of Ae are a measure for the 
strength of the NLO corrections in the effective potential. The 
units are [V] — [m„] = GeV. 



However, one can obtain the effective potential directly 
from the gap equation. Indeed, let us assume that the po- 
tential U ((Ta, 0a) of the Lagrange density of the bosonized 
theory is known, then < U{aa, <j>a) >= V{Ai). To explore 
the properties of the spontaneously broken theory, we re- 
strict ourselves to the part of the total potential, U{ai), 
involving only the fields which develop a nonzero vacuum 
expectation value, < ai >= Ai. Expanding U{ai) about 
the asymmetric ground state, we find 



(a. - A,) 



(67) 



It is clear that U{Ai) = V{Ai) 
functions of Ai 



and the derivatives are 



3.3 NLO contribution to the effective potential: SU{3) 
symmetric result 

The expectation value of the energy density in a state 
for which the scalar field has the expectation value Ai is 
given by the effective potential V{Ai). The effective po- 
tential is a direct way to study the ground state of the 
theory. If V{Ai) has several local minima, it is only the 
absolute minimum that corresponds to the true vacuum. 
A sensible approximation method to calculate V{Ai) is 
the semiclassical expansion ^7]. The first term in the ex- 
pansion of V is the classical potential. In the considered 
theory it contains the negative sum of all nonderivative 
terms in the bosonized Lagrange density which includes 
the one-loop quark diagrams and the leading order SPA 
result (|30|l . The second term contains semiclassical cor- 
rections from Ea. H46|l and the one-loop meson diagrams. 



dU 
da. 



dU 
dai 



dA, 



hi^j)- (68) 



It means, in particular, that we can consider Ea. (|68|l as a 
system of linear differential equations to extract the effec- 
tive potential V{Ai), if the dependence fi{Aj) is known. 

Further, Eq.jnH) tells us that fi{Aj) are determined 
by the tadpole term in a shifted potential energy, ?7((T^ -|- 
Ai), where we define a new quantum field with vanishing 
vacuum expectation value a^ = ai — Ai. In the case of 
SU{3) fiavor symmetry, we have fi = f where i = u,d,s 
and / is given by 



2 47r2' 

Kn-^uj{i - 3^2) 

4G2(l + 2tj)(l-4tj2)(i_t^2 



(69) 
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Therefore, the condition for the extremum dV{Au)/ dAu = 
coincides with the gap equation (|48|l . In Appendix ICl we 
obtain from Ea. H68|l in the SU{3) limit the effective po- 
tential 



Mml) 



-n 



4 "7 2 
-iln|(l-4cj2)(i_cj2)8| -a (70) 



The free constant C can be fixed by requiring 1^(0) = 0. 
In this expression w(m^) is defined according to Ea. ljl^lll 
and hu is the first of the two solutions of the stationary 
point equation given by 152() . Notice, that these solutions 
are complex when 4G^ < kAu- Hence the function V{Au) 
is not real as soon as the inequality is fulfilled. The most 
efficient way to go round this problem and to define the 
effective potential as a real function on the whole real axis 
is to treat as an independent variational parameter in- 
stead of Au in Ea. (|7()(l . In this approach, which actually 
corresponds more closely to the BCS theory of supercon- 
ductivity, one should consider equation (|126|l as the one 
which yields the function Z\,j(/i,j). One can check now that 
the extremum condition dV{hu)/dhu = is equivalent to 
the gap equation (|48() where the quark mass is expressed 
in terms of hu (16^ . Thus, the effective potential in the 
form of V{hu) provides for a direct way to determine the 
minimum of the vacuum energy irrespectively of the not 
well defined mapping A^ — > hu- 

There is a direct physical interpretation of Ea. H7U|) : 
classically, the system sits in a minimum of the poten- 
tial energy, Ucu determined by the first two terms, and 
its energy is the value of the potential at the minimum, 
Uci{Mu)- To get the first quantum correction, AU, to this 
picture, we add the third term (~ f^^^), and approximate 
the potential, V{Mu + Amu), near the classical minimum 
by a function Uci{Mu) + AU(Mu). 

In Figl^lwe show the effective potential calculated for 
K = -1800 GeV"^ G = 10 GeV"^ m = and yl = 
860 MeV depending on the strength of the fluctuations, 
indicated by the Euclidean cutoff yl^; . In absence of fluc- 
tuations the minimum occurs at m„ = Mmm = 340 MeV 
(outside the range of this figure). Increasing the effect of 
fluctuations, the minima Mmin appear at smaller values 
and the potential gets shallower. Simultaneously a barrier 
develops between V{Q) and V{Mmin)- At some critical 
value Ae the point V{Q) becomes the stable minimum 
and the trivial vacuum is restored^. This effect has a sim- 
ple explanation. In the neighbourhood of the trivial vac- 
uum where vriu is small the effective potential V^(m„) can 
be well described by the first terms of the series in powers 
of m,. 



3774 

4G 
0{ml 



NcGA^ 
27r2 



32G3 



Ae 
2tt 



(71) 



^ Away from the chiral limit the barier between the two 
vacua ceases fast to exist and the transition from one phase 
to the other occurs smoothly. 



The trivial vacuum always exists when 
A^^Gyl^ k2 [Ae 



1 - 



27r2 



32G3 V 27r 



> 0. 



(72) 



This inequality generalizes the well-known result for k — 
0. 

The local minima Mmin in the broken phase are the 
exact solutions of the full gap equation (|48|l . We find that 
at leading ?i-order, M^[^ = Af„ -I- Aniu, where Z\m„ 
is the correction 1)61(1 . follows within a few percent the 
pattern shown for Mmm in FigEI For instance, one has 
at Ae = 1.25 GeV, M^7„* = 214 MeV, and at Ae = 
1.31 GeV, M^[^ = 188 MeV. It is clear that the phase 
transition shown in FigOlis a non-perturbative effect. In- 
stead the perturbative result yields M^[^ smoothly 
with increasing Ae up to the value Ae — 1.6 GeV. 

Going to higher values of (not shown in Fig|3|) one 
can come to caustics, i.e. singularities in F(Z\„). From 
the logarithm in Eg. ((70)1 we obtain the values Z\„ where 
it happens. There are two singular points 



A. 



(1) 



ri2 r<2 
-12— , A(^) = -32— 



(73) 



For given values of couplings G and k we have = 
667 MeV and Z\lf ^ = 1.78 GeV. The indicated curves with 
Ae have as asymptote the vertical line crossing the 
mtj-axis at the point Au"^ . It is clear that for other pa- 
rameter choices the ordering A'^^^ < Mu < A'^\ or even 
Au'' < A^u^ < Mu are possible, where M„ is the classical 
minimum. In these cases a careful treatment of the caustic 
regions must be done. 



3.4 Gap equation at leading order: SU{2)i x U{1)y 
case, general properties 

We will now apply the same strategy to the case rhu = 
7^ TOs, which breaks the unitary SU{S) symmetry 
down to the SU{2)j x C/(l)y (isospin-hypercharge) sub- 
group. In full agreement with symmetry requirements it 
follows then that = nid ^ rris and hu = hd ^ hs- 
Thus, we have a system of two equations (from Eg. ((31(1 ) 
to determine the functions /i„ and h^. These equations can 
be easily solved in the limit G 



-Au\f^ (G^O). (74) 

K V As 



Obviously the result Eg. 1(53(1 follows from these expres- 
sions. For G 7^ the system is equivalent to the following 
one 



(75) 



Ghs + As 
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Fig. 4. Graphical solution of the cubic equation in 17511 . 
The left-hand (/i) and the right-hand (/2) sides of this equa- 
tion are plotted as a function of hu for Xs < 1. The region 
bounded by the dashed lines corresponds to the interval D < 
where the cubic equation has three real solutions. The lo- 
cal maximum A and the local minimum B have coordinates 
A:{h^ = -V^, f = A{xs)/G} and B : {h^ = ^/^, / = 
-Aixs)/G}. 



where we put Xg = kAs/ (AG)'^ in accordance with our no- 
tation in the Appendix 1X1 From the physics of instantons 
we know that the strength constant k < 0. It makes Xg 
negative. Hence the left-hand side of the first equation is 
not a monotonic function of /i„ (see Fig0J. The right-hand 
side is a positive constant, if Z\„ > and G > 0, what is 
usually assumed. Therefore, the equation has three dif- 
ferent real solutions, h^\Au, As), n = 1,2,3, in some 
interval of values for Au/G. The boundaries of the inter- 
val are given by the inequality ~A{xs) < Au < A{xs), 
where 



A{xs) 



32G2 



3/2 



(76) 



By means of the discriminant D of the cubic equation, 
D^Q^ + R'^, where 



16G\ Xs - 1 
~) 3 



(77) 



this region can be shortly identified by < 0. The qual- 
itative picture of the dependence /i„(Z\„) at a fixed value 
of As is shown in Fig[S| The three solutions in the region 
D < can be parametrized by the angle ip 



/iW =2v/^cos| , h^^^=2y 




(78) 



Fig. 5. The stationary trajectory hu{Au) at fixed As- Inside 
the interval \Au\ < ^(2:3) it is a multi- valued mapping. The 
monotonic curves: BC, CD, DE correspond to the three well- 
defined single- valued functions: ,n = 1, 2, 3. 



where 



COS if 



R 



simp ■ 




(79) 



The angle ip can always be converted to values of (p such 
that < ip < n. The boundaries (p = and ip> = it corre- 
spond to the value D ~ 0. When the argument (p increases 
from to TT, the solutions hu"^ , h^u"^ and h^u"^ run along the 
curves BG, DG and DE accordingly. These curves inter- 
sect the hu-SLxis at ip — n/2, where R — 0. One can show 
(2) 

that hu —Au/G at k ^ in full agreement with our 
previous result following from Ea. H52() . On the other side 

h^u^ leads to the result 1)74(1 in the limit G — > 0. Previously, 
studying the SU{Z) case, we have obtained both of these 
limits from one solution 1(521) . Now one has to use either 

solution h'^u^ or solution h'^u^ depending on the values of 
the parameters k, G, Au , and As which correspond to the 
minimum of the energy density. There is no chance to join 
the partial cases G — and k = in one solution, because 
they lead to systems of quadratic (in the first case) and 
linear (in the second case) equations without intersection 

of their roots. The stationary trajectory /li^^ is a positive 
definite function and thus the solution of the gap equation 
does not belong to this branch. 

Let us suppouse that the system of two equations which 
describe the vacuum state of the theory in the case of the 
SU{2)i X U{1)y symmetry at leading order has a solu- 
tion, i.e., the constituent quark masses Mui^ij and Ms(^i^, 
corresponding to a local extremum of the effective poten- 
tial V{'mu,ms), are known. We may assume that Mu(i) 
belongs to the region with the barely broken symmetry. 
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Mu(i) < fhu- As we already know, it is the most preferable 
pattern from the physical point of view. Then there is no 
any other solution -^^^(i) 7^ -^^m(i) for this set of param- 
eters G, K, A,Thu,Tris and already fixed value Mg(^iy This 
follows from the pure geometrical fact that the second or- 
der derivatives for curves and < QuQu > have op- 
posite signs in the considered region. Suppose further that 
there are other solutions M„(„-), Ms(„) with n — 2,3. This 
statement does not contradict our previous result corre- 
sponding to the case with the SU{3) symmetry, where 
we were able to find out three solutions for some sets 
of parameters. The functions M„(„) and M^^^n) can be 
understood as chiral expansions about the SU{3) sym- 
metric solutions. The coefficients of the chiral series are 
determined by the expectation values Qij{Mu) and their 
derivatives. It means that solutions in the case of SU{3) 
symmetry are related to the solutions for the more general 
SU{2)i X U{1)y case. Hence, there are only three sets of 
solutions (M„(„), Msj-n)) at maximum (for the considered 
region). The effective potential helps us to classify these 
critical points as will be discussed in Sec. 13.61 



3.5 NLO corrections: z, and Ami 

SU{2)i X U{1)y symmetry 



in the case of 



Our next task is to take into account quantum fluctuations 
and compute the corresponding corrections to the leading 
order result M^, Mg. For this purpose we need to find 2„ 
and Zs- Consider first the sum Aacthb which can be written 
as a 9 X 9 matrix in block diagonal form 



Aacbhb — 



i7i 

r?2 

!73 



(80) 



with 2 X 2, 3 X 3 and 4x4 blocks 



1 f-2{2K + hs) V2{K~hs) 



(81) 



The indicies r, s of the first matrix range over the subset 
r,s — 0, 8 of the set a,c = 0, 1, ... 8. In the matrix 
we assume that 71,111, = 1,2,3 and in the indices take 
values f,g = A, 5, 6, 7. 

Using this result one can solve the last two equations 
in H27|l . rewriting them in the form 



16G 



G 



V 



It 











16G 











hg'^^^-Sae (82) 



16G/ ac 

and find the functions hce ■ We obtain 



-6„ 



G{lT^s) ' 



, (1,2) ^ -Sfg 



(83) 



where are defined in the Appendix ^ For the 2x2 
matrix with indices 0, 8 we have 



3G(1± - 2ujI) 
3 =F (4a;„ - ujs) ±V2(a;„ - ujs) 

±V2(wm— Ws) 3 ± 2(2a;„ + o^s) ^ 



(84) 



In particular, if the terms lj„ and ojs are equal, these ex- 
pressions coincide with Eq.jS^. 



We now calculate the sum Aabchab^ where ha 



,(1). 



.(2) 



Using the properties of coefficients Aabc and solutions 



for /i^^'^'' obtained above, one can find that 



An brh, 



abc'^ab 



^(^88 - 2/100 
{hs8 + V2/108 - 



3%/3 
HiSco 



- 3/iii 
3hii - 



I- 4/l44)JcO 
2h44)Sc8 



HoS. 



20 ea- 



rn 



Again, from this equation the related formula (|57|l can be 
established by equating a;„ = uJs- In this special case we 
have H2 = 0. As a next step let us contract the result with 
functions h 



(1) 



A h - — 



HiiV2hil^ 



H2{V2h 



(1) 



"-88 /' 



(86) 
i — u,d 



Hi{V2h\^ - 2h\,'^) + H2{V2h\^ - 2h'i^), I = s 

These contributions can be evaluated explicitly. It leads 
to the final expressions for the semiclassical corrections Zi 
to the gap equation (|48|l . They are given by 



2(1 - 2ul 



8G2^ 



(l-2c.2)2- 
4c^2(3^^-2) 



il-2u;ir 
4 



1 - t^s 



1 



- OJf, 



(87) 

where /i = {l + ajg — 2ti>2). Each of these formulas has the 
same limiting value at a;^ = a;„ = a; which coincides with 
Eq.(|Sni). 

We apply this result to establish the ?i-order correction 
to the masses of constituent quarks. To this end we must 
use the general expressions obtained in Appendices IXI and 
IBl which for the considered case we rewrite in a way that 
stresses the quark content of the contributions 



Arm = -G ^ ZjMjkQ{i)k 

j,k=l 
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where i = u, s, the ?i-corrections Zj are written as a Une 
matrix Zj = (z„,Zs), and obviously Aruu = Arud- The 
2x2 matrix Ai jk and the column Q(i)k are defined as 



follows 



M 



1 + ujg 2uj 

UJ 1 



Q(i)k 



Qiu + Qid 



(89) 

Observe that detA^ = fi. The ?i-corrections to the quark 
masses must be calculated at the point (M„, Mg), being a 
solution of the gap equation at leading order. The expres- 
sion (|61|l is a straightforward consequence of the more gen- 
eral result (|88|l . Let us also note that formula H88|l clearly 
shows which part of the correction is determined by the 
strange component of the quark sea and which one by the 
non-strange contributions. 



3.6 Effective potential: the SU{2)i x U{1)y symmetry 

We can now generalize the result obtained in the Sec. 13.31 
to the case of SU{2)i x U{1)y symmetry. To find the ef- 
fective potential for this case one has to evaluate the line 
integral of the form 



(90) 



where the independent variables m„ and are linear 
combinations of the SU (3) singlet and octet components 

m„ = (v^Tio + m%)/^f^, iris = (v^wo — 2m8)/'\/3- We 
also know that 

fu[mu,ms) = — - - ^m„Jo(m„) + — , 
/^(m^,™^) = -y - -^msJo(m^) + y . (91) 

The functions fu, fs lie on the surface S{'mu, rris) defined 
by the stationary point equation l|75|l . and 7 is contained 
in S. There are some troubles caused by the singularities 
in Zu,Zs. The poles are located on curves which divide 
the surface S{mu,ms) on distinct parts i7„. The integral 
H9l)|l is well defined inside each of these regions i7„ . It is 
characterized by the property that the integral over an 
arc 7, which is contained in S^, depends only on its end 
points, i.e., the integral over any closed curve 7c contained 
in Sn is zero. This follows from the fact that the integrand 
is an exact differential. Indeed, one can simply check that 
the one- form 2fudmu + fsdrus is closed on Sn'- 



dfu dfs 



9m, dr 



= 0. 



(92) 



On the other hand, the open set Sn is diffeomorphic to 
E? and, by Poincare's lemma, the one-form is exact. 

Direct verification is relatively cumbersome and can 
be done along the lines of our calculation in the Appendix 
O as follows. Since each of the differentials 2/i„(im„ -t- 



hsdrus and 2zudmu + Zgdnis is closed, let us consider them 
separately. We begin by evaluating the first one-form 



{2hudmu + hsdnis) =d[Ghl + ^h^ + ^hlhs ] (93) 



where we have used Eas. (l31|) to extract drriu and drris. 

Noting that 3k^(/i^ + Q) = — /^(16G)^, one can obtain 
for the second one-form 



2zudmu + Zgdrris 



3k 
16 



{zs + 2uJuZu) dQ - 2GiiZudhu 



Let us note also that 



Zs + 2lOuZu 



(94) 



86-2 



1 



(1-2C.2)2-C.2 



Putting this expression in Ea. H94|) . we have after some 
algebra 

2zudmu + Zgdnis (96) 
duj] + 4(1 - 2ujl)dujl MujI MujI "I 



fl - 2C^2)2 



1 — U!-; 1~ Lot 



= -n-'d\n I [(1 - 2^2)2 _ ^2](i „ ^2)3(1 „ ^2)4| _ 

Finally, we obtain the effective potential 



V{A^,As) = UGhl + ^hl+'lhlh, 



-v{ml) - ]^v{ml) - C 



(97) 



\nm-2u:iY-ui]{l-u^iY{l-u^i) 



where v{mf) has been introduced in Ea. H131|l . The con- 
stant C depends on the initial point of the curve 7, and in 
the region S which includes the point Au = 0, Z\s = it 
can be fixed by requiring V{0, 0) — 0. This result coincides 
with Ea. fTlHl in the SU{3) limiting case. 

In FigElwe show the effective potential calculated as 
function of condensates without the fluctuations, ]7^^ = 
0, for the parameter set given in the caption. There are 
altogether nine critical points: four minima, four saddles 
and one maximum. Only one critical point is localized in 
the region of physical interest, the minimum for uiu, rus > 
(or hu,hs < 0). In the chiral limit and otherwise the 
same parameters one has three critical points of interest, 
the maximum at the origin, the saddle at ruu = 0, nis > 
{hu — O^hs < 0), and the minimum. This distribution 
and behavior of critical points is common to a large set of 
parameters. 

The behavior of V{Au, As) in terms of the strength of 
the fluctuation term ~ is qualitatively the same as 
for the SU{3) case: fluctuations tend to restore the trivial 
vacuum in the region prior to the singularities. We also 
see from Ea. H97(l that now the picture of singularities is 
more elaborated. Nevertheless attractive wells still develop 
between them. 
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Fig. 6. Classical effective potential l/[GeV], Eg. 19711 with 
— 0, as function of condensates [GeV""^] hu = 2 < uu >, 
h, = 2 < ss > ioi A = 860 MeV, G = 14 GeV'^ k = 
-1000 GeV-^ m„ = 6 MeV,m, = 150 MeV. 



These results are in agreement with the following topo- 
logical consideration. The effective potential F is a smooth 
function defined on the space of paths S^h^, hs) diffeomor- 
phic to (before the onset of caustics). The Euler char- 
acteristic of the surface S, xi'S) = 1, can be expressed, by 
Morse's theorem, through the number of non-degenerate 
critical points of the function V 



X{S) ^Co-Ci+C2 



(98) 



where Cq is a number of critical points with index (min- 
ima) , Ci is a number of critical points with index 1 (saddle- 
points), and C2 is a number of critical points with index 
2 (maxima). 



4 The ground state in the 1/Nc expansion 

We have considered till now the semiclassical approach 
to estimate the integral Z[a,4);A] in Eg. ljiyi) . However, 
rather than using h as the parameter of asymptotic ex- 
pansion, we could also have used 1/Nc to estimate it. In 
this case the stationary phase equations involve terms 
with different orders oil /Nc and must be solved perturba- 
tively. For this purpose let us represent (|21|l in a complex 
form 



GUa + Wa 



Sac 

32^' 



.bcUlUl 



0. 



(99) 



Indeed the first two terms are of order (l/A^^)", since G ~ 
1/iVc, U ^ Nc, and W ~ iV°, while the last one is of order 
l/Nc- Casting the solutions U^^ as a series in 1/Nc up to 
and including the terms of order 1/Nc, we have 



Tja 

^st 



3k 
32G2 



Aa 



be 



wlwl 



0{1/Nl) 



(100) 



It yields for Cr{r^\/) 

+ 0{1/Nc). 



(4G) 



3 (det W + det W"^) 

(101) 



The contribution from the auxiliary fields can be ob- 
tained directly from Ea. H47|l by expanding our solutions 
hab and he in a series in l/Nc- One can already conclude 
from that expression, without any calculations, that the 
term ~ fi~^ is at most of order ~ l/Nc- Therefore it is 
beyond the accuracy of the considered approximation and 
can be neglected. Actually, it follows from Ea. (|27|l . that 



the difference hab — ha 



(1) 



''ab 



has order ~ iV°, addition- 



ally suppressing this contribution, i.e. for the correction to 
the result pOlfl we have only a term starting from 1/Nc- 
order 



~8(2G)4 

o{i/n!). 



ti{WW^) 



3k 

JagY' 



(det W - 



det W'') 
(102) 



It means, in particular, that one can neglect this type of 
quantum fiuctuations in the discussion of the gap equa- 
tion up to and including the terms of l/Nc order in the 
meson Lagrangian. These corrections cannot influence sig- 
nificantly the dynamical symmetry breaking phenomena 
in the model and only the quantum effect of mesons (the 
one-loop contributions of a and (j) fields) together with the 
leading order contribution from the 't Hooft determinant 
(see Ea. (|101|) ') are relevant at A^"-order here. 

If the model allows us to utilize the l/Nc expansion, 
the vacuum state is defined by the gap equation obtained 
from Ea. (|48l) in the large Nc limit, or, equivalently, on the 
basis of Lagrangian H101|) . We have 



A,, 



32G2 



N G 

UjkAjAk = -^-^■miJoimi). (103) 



At leading order in an 1 /Nc expansion we have the stan- 
dard gap equation ^iP' Ai = NcGmiJQ{mi). The terms 
arising on the next step already include the quantum cor- 
rection from the 't Hooft determinant. It is not difficult 
to obtain the corresponding contribution of order 1 /Nc to 
the mass of constituent quarks 



Am,: 



M, + Ar 



nAjAk 



(4G)2(l-/,0 



(i ^j^k) (104) 



where Ami is calculated at the point ~ Mi, which is 
the solution of the gap equation H103|) at leading order, 
and li is given by Ea. Ht)3l) . One can show that Ami > 0, 
thus increasing the effect of the dynamical chiral sym- 
metry breaking. This is an immediate consequence of the 
formula 

1 - ^^(*^') = ^ + ^MfMM^) > 0. (105) 



16 



A. A. Osipov and B. Hiller: Path integral bosonization of the 't Hooft determinant: quasiclassical corrections 



Let us stress that the first equality is fulfilled only at the 
point Mi. 

We did not clarify yet the counting rule for the current 
quark masses rhi, assuming that they are counted as the 
constituent quark masses. Actually, these masses are small 
and following the standard rules of ChPT one should con- 
sider rhi ^ i/Nc- In this case in the large Nc limit the 
model possesses the C/(3) x U{3) symmetry and the gap 
equation at leading order, 27r^ = NcGJo{mi), leads to a 
solution with equal masses = Md = Mg = M . The 
l/Nc correction includes the K-dependent term and the 
term depending on the current quark masses 



Ami — 



kNP \ 1 
16G2 ) M^Ji{AP) 



(106) 



Returning back to Ea. (|101(l one can conclude that the 
large- iVc limit corresponds to the picture which is not af- 
fected by six-quark fluctuations. This can be realized also 
directly from Lagrangian (|30|l . Indeed, the couplings of 
meson vertices are determined here through the functions 
ha given in the case of SU (3) flavor symmetry by h'^u^ in 
Eg. 1)52(1 . The large- A^c hmit forces a series expansion for 
hu"^ with a small parameter 



\k\Au 
,(1) - 



1 



< 1 



-Au/G^ 



(107) 



which does not 



and the leading term hi^ 
depend on k. 

This observation leads us to the second important con- 
clusion. It is easy to see that the parameter e is an inter- 
nal model parameter and the series expansion in e closely 
corresponds to the 1/A''c expansion of the model. The ex- 
istence of this small parameter allows us to consider the 
1/Nc series as a perfect approximation for the system with 
small vacuum six-quark fluctuations 



|k| « 



4G^ 
Au 



A^ = Au{G,A). 



(108) 



What to do if e is not too small? A large value for e can 
simply distabilize the 1/Nc series, implying large 1/A^c cor- 
rections. It has been observed recently [TH! that the abun- 
dance of strange quark-antiquark pairs in the vacuum can 
lead to nonnegligible vacuum correlations between strange 
and non-strange quark pairs. If this happens one can try to 
understand the behaviour of the quark system on the basis 
of the ?i-expansion. This approximation can be considered 
as the limit of large six-quark fluctuations. Although QCD 
does not contain an obvious parameter which could allow 
one to describe this limit, the model under consideration, 
as one can see, for instance, from Ea. lfT^ . suggests this 
dimensionless parameter: 



32G3 



< 1. 



(109) 



On the basis of this inequality one can conclude that val- 
ues of K, corresponding to large six-quark fluctuations are 
determined by the condition 



where we used letters I, s 
k\ for large and small six-quark 



One can see that ~ y/Nc 
to mark possible values of 
fluctuations correspondingly. These two regimes lead to 
the different patterns of chiral symmetry breaking. In the 
first case the mass of 77' meson goes to zero when Nc — > 00. 
In the second case the ground state does not have the r]' 
Goldstone boson even at leading order. 



< 2v/(2G)3|2. 



(110) 



5 Concluding remarks 

The purpose of this paper has been to use the path in- 
tegral approach to study the vacuum state and collective 
exitations of the 't Hooft six quark interaction. We started 
from the bosonization procedure following the technique 
described in the papers 6 ,8 . The leading order stationary 
phase approximation made in the path integral leads to 
the same result as obtained by different methods, based 
either on the Hartree - Fock approximation [I], or on the 
standard mean field approach |16) . The stationary trajec- 
tories are solutions for the system of the stationary phase 
equations and we find them in analytical form for the two 
cases corresponding to SU{3) and SU{2)i x U{\)y fla- 
vor symmetries. The exact knowledge of the stationary 
path is a nesessary step to obtain the effective bosonic 
Lagrangian. We give a detailed analytical solution to this 
problem. 

As the next step in evaluating the functional integral 
we have considered the semiclassical corrections which 
stem from the Gaussian integration. We have found and 
analyzed the corresponding contributions to the effective 
potential V , masses of constituent quarks and quark con- 
dencates. The most interesting conclusions are the follow- 
ing: 

(a) We have found that already the classical effective 
potential V{Au) for the case in which chiral symmetry 
SUl{3) X SUr{3) is broken down to the SU{3) subgroup, 
has a metastable vacuum state, although the values of 
parameters G, K,rhi, A corresponding to this pattern are 
quite unnatural from the physical point of view: the cou- 
plings G and A must be small to fulfill the inequality 
NcA^G < 27r^, which is known in the NJL model without 
the 't Hooft interaction as a condition for the trivial vac- 
uum; the coupling k must be several times bigger (in abso- 
lute value) of the value known from the instanton picture. 
Besides, the window of parameters for the existence of 
metastable vacua is quite small. We show then that semi- 
classical corrections, starting from some increasing critical 
value of the strength Ae, transform any classical potential 
with a single spontaneously broken vacuum to the semi- 
classical potential with a single trivial vacuum. Close to 
the chiral limit this transition goes through a smooth se- 
quence of potentials with two minima. There arc other 
known cases of effective chiral Lagrangians which confirm 
the picture with several vacua [TUll^ . 

(b) If the symmetry is broken up to the SU{2)i x 
U{1)y subgroup the smooth classical effective potential 
V{Au, As) defined on the space of stationary trajectories 
S may have several non-degenerate critical points. It is 
known that the properties of the critical points are related 
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to the topology of the surface S. We used this geometri- 
cal aspect of the problem to draw conclusions about an 
eventual more elaborate structure of the hadronic vacuum 
already at leading order in h. We find for some parameter 
sets the existence of a minimum, a maximum and a saddle 
point. For vanishing current quark masses, for example, 
the minimum corresponds to the spontaneous breakdown 
of chiral symmetry, the maximum is at the origin, and 
the saddle at m„ ^ and finite. Similar as in the 
SU (3) case, the inclusion of fluctuations tends to destroy 
the spontaneous broken phase and to restore the trivial 
vacuum. 

Our work raises some issues which can be addressed 
and used in further calculations: 

(1) The Gaussian approximation leads to singularities 
in the effective potential. To study this problem which 
is known as caustics in the path integral, it is necessary 
to go beyond this approximation and take into account 
quantum fluctuations of higher order than the quadratic 
ones. To the level of accuracy of the WKB approximation 
the effective potential is elsewhere well-defined and has 
stable minima between the singularities. It is interesting 
to trace the fate of the these minima going beyond the 
WKB approximation, since in this work we have analyzed 
the effective potential mainly at a safe distance from the 
first caustic. 

(2) Our expressions for the mass corrections (|117|l have 
a general form and can be used, for instance, to include the 
one-loop effect of mesonic fields. One can use it as well to 
find the relative strength of strange and non-strange quark 
pairs in this contribution. 

(3) We have chosen h and independently 1/Nc as two 
possible parameters for the systematic expansion of the 
effective action. As we saw the l/Nc expansion is a much 
more restrictive procedure. They are of interest for the 
study of the lowest lying scalar and pseudoscalar meson 
spectrum. There is a qualitative understanding of this 
spectrum at phenomenological level (see, for instance, pa- 
pers 12]). A more elaborate study might lead to the neces- 
sity of including either additional many-quark vertices or 
taking into account systematically quantum corrections. 
Our results might be helpful in approaching both of the 
indicated developments. 
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following from Ea. H48|l 



Edhu . 



-Arrir. 



Id 
G' 



A 

'g 



i—u,d,s 

E 



dhd 
dA, 



i—u,a,s 



Ami 



At 



Zd 



(111) 



where the functions Ii{mf) are given by Ea. (|63|l . Both the 
partial derivatives and integrals liimf) must be calculated 
for mi = Mi. From the first system of equations in (|50|l we 
may express derivatives dhi / dAj in terms of the functions 
hi. To symplify the work it is convenient to change the 
notation and rewrite H31|l in the form 



UJu + Xu+ i^di^s = 
i^d + Xd+ UJuUJs = 
UJs+Xs+ LOuiVd = 



nhi 
T6G 



nAi 
(4Gj2 



(112) 

Straightforward algebra on the basis of these equations 
gives 



dxi 



1 duj 
-(uJiUjj -uj-k)- — - 



A 



dxi 
(113) 

Here A = 1 — uj-^ — uj-^ — uj^ + 2L0u^d^s ■ We assume that 
indices J, j, fc range over the set {u, d, s} in such a way that 
* 7^ J 7^ ^ J 3' sum over repeated index is implied only 
if the symbol of the sum is explicitly written. 

The main determinant of the system (|111() is equal to 



P = - 



D 
AG3 



D = \- 



i—u,d,s 



h + (1 - u;l)ldls 



+ (1 - iuj)IJu + (1 - u;^s)luld - Al^hls 



(114) 



with li given by Ea. H63(l . The other related determinants 
are written in the compact form 



V,, = 



where 



AG2 



B: 



{ziBjk + ZjBjik + ZkBkij) , 



(115) 



{l-Lof)Ij-il-u;])I, + AIJj, 



Bijk = (t^iWj - uJk)Ik + ^/c- 



(116) 



Hence the ?i-correction to the mean field value of the con- 
stituent quark mass is given by 



A The semiclassical corrections to the 
constituent quark masses 

To derive the explicit formula for the semiclassical next 
to the leading order corrections Ami to the constituent 
quark masses Mi one has to solve the system of equations. 



D G 

Am^ = ^ = {ziB.jk + ZjBjik + ZkBkij) (117) 

where i ^ j ^ k. This formula gives us the most general 
expression which has to be specified by the explicit form 
for Zi. 

Let us also write out the two partial cases for this 
result. If Lu-u ~ LUd = '-^ and ^ Id = I, which happens 
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when the group of flavor symmetry is broken according to 
the pattern SU{3) SU{2)f x ^7(1)/, one can find 

A = {l-ujs)fi, /!= -2^2), 

D=[l-{1- CJ,)/][1 -/,-(!+ LU,)I + ^llls], 

AG^Vs = [l~{l-Us)I][zAl~ ^^I) + 2zuLo]. (118) 

This determines the coefficients in Ami, and we obtain 
the following form of mass corrections 



Z\m„ = G 

Anis = G 



Zuflls - Zu{l + LPs) - ZsOJ 

l-Is-{l+ uJs)I + nils 

Zsjfll - 1) - 2ZuLU 
1 - /, - (1 + UJs)I + fills 



(119) 



The second partial case for the formula H117|) corresponds 
to the SU{3) flavor symmetry. It is clear that now we have 

A = {l + 2uj){l-uj)\ 

D= [l-(l + 2cj)/][l-(l-tj)/]2, 

AG^Vu = z„(l + 2w)[l - (1 - Lj)I]^. (120) 

Then (|117|l yields the following result for the mass correc- 
tion 

Zu{l + 2uj) 



Arriu — —G 



1 - (1 + 2uj)l 



121) 



B Particle expectation values 

The knowledge of the constituent quark mass Mi gained 
from the equations (|50|l . combined together with the Feyn- 
man - Hellmann theorem have been used for finding 
the expectation values of the scalar quark densities, gj^j, 
in the physical quark state \Qi > (see paper ;4 ). The 
matrix element < QilqjqjlQi > describes the mixing of 
quarks of flavor j into the wavefunction of constituent 
quarks of flavor i. One can determine these particle ex- 
pectation values by calculating the partial derivatives 



< 



\QjQj\Q^ >-- 



drh-j 



(122) 



The functions Qij have been calculated in Q for the case 
m„ = md- There are a number of physical problems in 
which these matrix elements are usefull. We have found 
the presence of them in the expression for the semiclassical 
corrections to the mean field quark masses Mi. 

Both the mass corrections Anii and the particle expec- 
tation values Qij are solutions of a similar system of equa- 
tions and can be derived on an equal footing. Indeed, Qij 
are the solutions of the equations obtained from Ea. H5()|l 
by differentiation with respect to rfij . These equations dif- 
fer from the system ()lllll only up to the replacements of 
variables Ami — > Qij and Zi dhi/dAj. Therefore we 
have 



Q^ 



V, 



V 



(123) 



With this way of writing the determinants 2?^^) we wish 
to stress that one can simply obtain them from "Di in 
Ea. (|115() through the above replacements Zi dui/ {Gdxj). 
The main determinant T) is not changed. By use of these 
formulas we are led to the explicit expressions 



Q^^^-[{l-I,){l-h)^U;^J,h] 



Qi] = J^I^B^Jk, i¥'j^k (no sum) (124) 

which correspond to the most general case m„ ^ 7^ 
rhs- The notation have been explained in Appendix 



C Effective potential 

The models which are considered here lead in the most 
general case to the potential V{Au, A^, Ag) = V{Ai). 
This function is a solution of Eq. (|68l) . One can reconstruct 
V{Ai) by integrating the one- form 



av = — amu + -7: dma 



dm,, 



dmd 



dV 
dm. 



dms = ^ fidmi. 

i—u,d,s 

(125) 

Thus, the derivation of the effective potential in the /„ = 
fd — fs case is simply a question of representing the gap 
equation f{mu) = in the form of an exact differen- 
tial 3/(TO„)(im„ — dV{mu), where m„ — ■y/2/ 3mo . One 
should also take into account the constraint 



GK 

which implies that 



(126) 



dmu = - ( G + -hu I dhu = 



16G2 



(1 + 2Lu)dLu. (127) 



Using this result one can obtain for the first term in 1)69(1 , 
for instance, 



hudmu = hu{G+ 



dh,, 



d 



G, 
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(128) 



A similar calculation with the third term in 16911 gives 
kQ^^ a;(l — 3uj'^)dmu 

2G2 (l + 2cj)(l-4t^2)(l_^2) 
(1 - 30)2)^^2 



-412" 



{1 - Alu^){1 - uj^) 



1 



= ^i7-Mln |(1 - 4uj^){l - uj^f\ 
3 



(129) 



To conclude the procedure we must then add to the effec- 
tive potential V{Ay^) the corresponding contribution from 
the second term 



27r2 



muJQiml)dmu = -dv{m^) 



(130) 
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where we have defined 



( 2. 



(131) 



All this amounts to calculate the effective potential in the 
form given by Ea. H7U|) . 
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